clc;
clear;
f=inline('2/sqrt(p1)*exp(-x.^2)','x')


f=@(x)2/sqrt(pi)*exp(-x.^2);
y=integral(f,0,1.5)

syms x,
I0=vpa(int(2/sqrt(pi)*exp(-x^2),0,1.5))